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The Radial Velocity Spectrograph (RVS) on board of Gaia needs to be calibrated using stable reference stars known in advance. The 
catalogue presented here was built for that purpose. It includes 1420 radial velocity standard star candidates selected on strict criteria 
to fulfil the Gaia-RVS requirements. A large programme of ground-based observations has been underway since 2006 to monitor these 
stars and verify their stability, which has to be better than 300 m s _I over several years. The observations were done on the echelle 
spectrographs ELODIE and SOPHIE on the 1.93-m telescope at Observatoire de Haute-Provence (OHP), NARVAL on the Telescope 
Bernard Lyot at Observatoire du Pic du Midi and CORALIE on the Euler-Swiss Telescope at La Silla. Data from the OHP and Geneva 
Observatory archives have also been retrieved as have HARPS spectra from the ESO archive. We provide a mean radial velocity in 
the SOPHIE scale for each star, derived from the combination of velocities measured with those instruments, after having carefully 
estimated their differences in zero points. In total, 10214 radial velocity measurements have been obtained for the 1420 stars. With 
a mean time baseline of 6.35 years, 92.9% of the candidates fulfil a target stability criterion of 300 m s . Three hundred forty-three 
stars are found to be constant at the level of 100 m s _1 over 10 years. Comparisons with earlier catalogues show excellent agreement 
for FGK stars, with zero-point differences lower than 100 m s _I and a remarkably low RMS scatter of 33 m s _I in one case, suggesting 
that the precision of the catalogue presented here is better than this value. This catalogue will likely be useful for other large-scale 
C3 ' spectroscopic surveys, such as APOGEE, Gaia-ESO, HERMES, and LAMOST. 
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1 . Introduction about five million stars down to V - 13 and iron and a-elements 

_ . , „ , „ _ abundances for about two million stars down to V - 12. Such 

Gaia, the next astrophysical mission of the European Space a , spectroscopic survey will have a tremendous impact on 

, Agency (ESA), will be launched in 2013U It will survey the sdence caseSi such afj ^ chemist and dynamics of the 

. entire sky during five years, and measure the positions, proper MOk w ^ detection and characterization of multi le 

O ■ motions and parallaxes of some 10? stars down to a visual mag- tems> and yariable ^ More of ^ expected 

rn mtude of 20, with a n expected paralla x accuracy of 10 to 25 ^as ieM from ^ Rys m described in Wilkinson et aL (200 5). 
t 7j ■ at 15 magnitude (deBrunne, 2012). Stars classified as single 
J> ' will have their effective temperature, surface gravity, and metal- 

licity estimated, as well as interstellar extinction. xhe RVS has a reso i v i ng powe r of 1 1500 and covers the 

The Radial Velocity Spectrometer (RVS) on board Gaia will spectral range 847 _ 874 nm ^ which includes the near-infrared 

; provide radial velocities (RV) of about 150 million stars down Can uiplet ^ many Unes of iron and a _ e l e ments, and Paschen 

to 17 l magnitude with precisions ranging from 15kms at the lines in hot stars xhis spec tral range has already proven to 

faint end to 1 km s' 1 or better for Gand K^satmagmtudes be very well su ited to RV analysis with the RAVE survey 

brighter than V - 13.5 (|Katz et aU l200j|Katil2009[). The RV, dSteinmetz et all l2006h . The RVS is an integral field spectro- 

combined with astrometry, will give access to the six dimensions graph with no entrance s i it and no on -board wavelength calibra- 

of the phase space useful for kinematical studies. The RVS will tion source As a conse quence, the RVS wavelength dispersion 

also provide rotational velocities and atmospheric parameters for law win be derived from the RVS obse rvations. The ins t rumen t 

I _ " 7~ T~- 7~ ■ i ■ i ^ • i — *n a i — i /t~\ i will be self calibrated. As explained in Jasniewicz et al. (201 1), 

Based on data obtained within the Gaia DPAC (Data Processing and , ... , . , ' : ' , 

Analysis Consortium) and coordinated by the GBOG (Ground-Based ground-based radial velocity standards are mandatory for cal- 

Observations for Gaia) working group, at various telescopes; see ab- Crating the wavelength zero point, which otherwise would be 

stract mathematically degenerated with the radial velocity scale (i.e. 

** Tables are only available in electronic form at the CDS via anony- a shift in the wavelength scale can be compensated for by a 

mous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via http://cdsarc.u- shift in the RV scale). The Gaia Data Processing and Analysis 

strasbg.fr/viz-bin/qcat ?J/A+ A/552/ A64 Consortium (DPAC), and more specifical l y its Development 

http://www.rssd.esa.int/index.php?project=GAIA&page=index Unit 'Spectrocopic Processing' (Katz et all 1201 ll) . is responsi- 
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ble for establishing the requirements for this external calibration 
and for the acquisition of the necessary auxiliary data 

In principle asteroids are excellent sources for calibrating 
the zero point because their radial velocity can be derived 
from celestial mechanics with uncertainties lower than lms 4 . 
However, they cannot be the main reference sources for Gaia 
because they are not numerous enough in the appropriate mag- 
nitude range and their sky distribution is limited to the Ecliptic 
Plane. 

When we started to search for suitable RV standard stars six 
years ago, there was no existing dataset fulfilling the Gaia RVS 
requirements in terms of sky coverage, magnitude range, number 
of targets, and stability. The only official list of standard stars 
was the one compiled by the IAU Commission 30, which did not 
turn out to be suitable for our purpose. We thus built a stellar grid 
consisting of 1420 FGKM stars in the magnitude range 6 < V < 
11, with a homogeneous distribution on the sky. The s e lectio n 
process of that sample is fully explained in ICrifo et al.1 (|2010). 
We briefly summarize it in Section [2] 

The RVS calibration requires standard stars with much bet- 
ter RV stability than the 1 km s _1 accuracy expected from then- 
future measurements with RVS, with no drift until the end of the 
mission (2018). The value of 300 m s _1 is adopted as the stabil- 
ity level to be checked. To qualify as a reference star, each candi- 
date has to be observed at least twice before launch and another 
time during the mission to verify its long-term stability. The RV 
measurements available to date are described in Section [3] and 
their combination into a homogeneous scale in Section [4] We 
compare this new catalogue to other studies in Section [5] The 
stablest stars are presented in Sect. [6]as are considerations about 
variable stars. Future work is described in Sect. [7] 

2. Radial velocity standard-star candidates 

The selection of stars suitable for building the reference grid 
for R VS has been caref ully considered and several criteria de- 
fined (ICrifo et all 12010). The list of candidates was established 
from three catalogues : "Radial velocities of 889 late-type stars" 
dNidever et all 120021) "Radi al velocities for 6691 K and M 
giants" dFamaevetafl, 120051) . a nd "The Geneva-Copen hagen 
Survey of Solar neighbourhood" (Nordstro m et all 12004). com- 
plemented with IAU standards dUdrv et all 1 19991) . We selected 
the stars having the best observational history in terms of con- 
sistency of radial velocity measurements over several years. In 
that way we could focus on stars more likely to be stable over 
time. Moreover, since all the candidates are in the HIPPARCOS 
Catalogue we were able to check their properties (photometry, 
spectral type, variability, multiplicity, etc.) in a homegeneous 
way. Another important criterion was that each candidate had to 
have no close neighbour within a circular region of 80"-radius. 
This criterion was set to ensure that their future RVS spectra will 
not be contaminated by another overlapping spectrum, 80" be- 
ing the length of an RVS spectrum projected on the sky. This 
criterion was verified thanks to the USNO-B1 catalogue. These 
different criteria led to a preliminary list of 1420 candidates hav- 
ing a high probability of being RV-stable and well suited to be- 
coming reference stars for the RVS. 

All of the 1420 candidates have at least three previous 
RV measurements in the above-mentioned catalogues; however, 
since they have been obtained during observational runs span- 
ning irregular time baselines, their stability has not yet been en- 
sured, particularly not until the end of the Gaia mission. The 
main cause of RV variability at the level of 300 m s _I and above 
is binarity, but the risk of having binary stars with significant 



amplitude in our sample has been reduced because of all these 
selection criteria and a careful analysis of the RV measurements 
available for them so far. The HR diagram of the 1 420 candi- 
dates is shown in Fig. Q~l as originally published in ICrifo et all 
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Fig. 1. HR diagram of the 1420 selected stars, with a colour code 
corresponding to their origin. 



3. Radial velocity measurements 

The RV standards to be used for the RVS calibrations should not 
exhibit variations larger than 300 m s _1 during the five years of 
the Gaia observations, between 2013 and 2018. Since the stan- 
dards need to be known before launch, we extrapolated this crite- 
rion to the current period in order to perform a pre-qualification 
of the candidates according to their stability during the time 
span of our ground-based observations, between 1995 and 2012. 
Considering that the candidates already have a good observa- 
tional history, our strategy of pre-qualifying them is (i) to obtain 
at least two new measurements per star before the Gaia launch, 
separated by more than one year, and (ii) to verify that the vari- 
ation in these new measurements does not exceed 300 m s _1 . 
Then supplementary ground-based observations will have to be 
done during the mission to verify the long-term stability of the 
reference grid. 

Since the candidates are distributed all over the sky, ob- 
servations in both the northern and southern hemispheres need 
to be performed with different instruments. An observing plan 
has been prepared under the umbrella of the GBOG (Ground 
Based Observations for Gaia) working group in charge of ac- 
quiring auxiliary data for the Gaia processing. Eighty nights 
have been obtained over six years for that programme on four 
high-resolution spectrographs: ELODIE then SOPHIE mounted on 
the 1.93-m telescope at Observatoire de Haute-Provence (OHP), 
NARVAL on the Telescope Bernard Lyot at Observatoire du Pic 
du Midi, and CORALIE on the Euler-Swiss telescope at La Silla. 

RV measurements were also retri eved from the publi c 
archives of SOPHIE and ELODIE at OHP dMoultaka et all 12004). 
and of HARPS at the ESO Advanced Data Products database. 
Many measurements are available in these archives, in partic- 
ular due to the various exoplanet-hunting programmes. The full 
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ELODIE archive was not public when the observations started in 
2006, and the volume of the other archives increased dramati- 
cally during the course of our programme. The recent query of 
the archives allowed us to retrieve several thousands of relevant 
measurements. There were many time series available for some 
stars. For time series in one single night, only the observation 
with the highest signal-to-noise ratio was considered. Additional 
data from the private CORALIE archive have also been provided 
by the Geneva Observatory. 

SOPHIE, CORALIE, ELODIE and HARPS have respectively re- 
solving powers of R = A./AA,- 75 000, 50 000, 42 000, and 
120 000, respectively, and cover the visible range. They all have 
similar automatic on-line data reduction software to derive the 
barycentri c RV by cross-corre lation of the spectra with a numer- 
ical mask (Bar anne et al.LlT996h . The spectral type of the numer- 
ical mask is chosen to be the closest to the observed star. It can 
be G2, K5, M4, and M5 for SOPHIE, G2, K5 and M2 for HARPS 
and CORALIE, F0 and K0 for ELODIE. 

NARVAL has a resolving power of R ~ 78 000. Since this 
intrument has not been built for RV measurements but for po- 
larimetry, the on-line reduction software does not include the 
RV determination. We measured it by cross-correlating the ob- 
served spectra with the SOPHIE G2 mask. Our main interest in 
observing with NARVAL is that its spectral coverage includes the 
RVS range (847-874 nm), allowing us to investigate possible 
systematic differences occurring when measuring the RV in the 
full visible range and in that narrow NIR spectral range. It is, 
however, beyond the scope of this paper to enter into these con- 
siderations. 

Errors on radial velocities come from the photon noise, the 
wavelength calibration and some possible instrumental or other 
effects, which are impossible to quantify. It is, however, impor- 
tant to estimate errors, noted e in the following, for the weight- 
ing scheme adopted to combine measurements taken with dif- 
ferent instruments and under various observing conditions. The 
photon-noise uncertainty is used, deduced from the width and 
contrast of the cross-correlation function (here after CCF) and 
from the signal-to-noise r atio of the spectrum (Ba ranne et al.L 
1996: iBouchv et all 120051) . and provided by the data reduction 
software. External systematic errors are then quadratically added 
to the photon-noise uncertainty: 0.8 m s _1 f or HARPS (F. Bouchy , 
private communication), 4 ms 4 for SOPHIE ( Boisse et alll2012l) . 
5 m s _1 for CORA LIE (D. Naef, private communication), and 15 
m s _I for ELODIE (iBaranne et al.. 1996). These systematic errors 
have been estimated for the most precise observing mode of the 
instruments where the stellar spectrum is recorded simultane- 
ously with a thorium-argon calibration. But not all the HARPS, 
SOPHIE, and ELODIE spectra in the sample presented here have 
been obtained in that mode, and the corresponding RV mea- 
surements have larger errors. This is taken into account by dou- 
bling the systematic error listed above for exposures without the 
simultaneous calibration. It is worth recalling that e estimated 
that way does not represent the total RV uncertainty, but is only 
meant to combine relative measurements with proper weights. 

Table Q] summarizes the status of the observations, provid- 
ing the number of spectroscopic measurements obtained so far 
with each instrument (Nrv), the corresponding number of stars 
(N„), the average errors of individual RV measurements e, and 
the date range of the observations. In total there are 10214 RV 
measurements. 

Figure [2] represents the distribution of the 1420 candidate 
stars on the celestial sphere, with a colour code indicating the 
number of measurements obtained for each star. As expected, at 
least two measurements are available for the vast majority of the 



Table 1. Summary of RV measurements available as of 
November 2012, average of individual RV errors, and epoch 
range of the measurements. 





Nrv 


N, 


6 


date 








m s -1 


range 


SOPHIE 


2945 


729 


6.7 


2006 - 2012 


CORALIE 


2470 


775 


12.9 


1999 - 2012 


NARVAL 


213 


157 


17.6 


2007 - 2012 


ELODIE 


3931 


372 


17.4 


1995 - 2006 


HARPS 


655 


113 


2.2 


2003 - 2009 



stars. One measurement is missing for only a few stars in the 
southern hemisphere. 

Figure [3] shows the histogram of the time span between the 
first and last measurements of each of the 1420 candidates. This 
is 6.35 years on average (median 4.37 years), up to 17 years 
thanks to the ELODIE archive. 
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Fig. 3. Distribution of time baselines of observations for the 
1420 candidates. 



4. Combining RVs from different instruments 

An important step in the process of determining whether a star 
has a stable radial velocity or not is to combine the measure- 
ments obtained with the different spectrographs. This requires 
taking the individual zero point of each instrument into ac- 
count. The zero point is related to the instrument itself, its spec- 
tral range, and resolution, and to the calibration procedure and 
the method used to derive an RV. In our case the same algo- 
rithm of cross-correlation with a numerical mask has been used 
that should minimize the offsets between instruments. However, 
since different masks have sometimes been used, a dependency 
of the offsets on the colour is expected. We adopted the SOPHIE 
frame as the radial velocity reference scale for all spectrographs, 
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Fig. 2. Distribution of the 1420 stars on the celestial sphere in equatorial coordinates. The colour code indicates the number of 
measurements obtained for each star. A dashed line indicates the projection of the Ecliptic plane, a dotted line that of the Galactic 
plane. 



and we made certain to have a sufficient number of stars ob- 
served in common with SOPHIE and with the other instruments 
to derive the offsets. 

To measure the offset between SOPHIE and the other instru- 
ments, we first computed the weighted mean and standard devi- 
ation of the RV measurements per instrument and per star. The 
weight of each individual measurement RV, is defined as 1/ef, 
e, being the RV error mentioned in Sect. [3] 

Figure |4] shows the RV offsets between the various instru- 
ments versus the B — V colour index. The error bars, when 
available, are the quadratic sum of the two standard deviations 
obtained for a given star observed several times with the two 
considered instruments. A large error bar is an indication that 
the star is probably variable. The offsets are constant in gen- 
eral, except for ELODIE where a dependency on colour is clearly 
seen. The scatter is also much higher when ELODIE is involved. 
The scatter reflects the RV variations, expected to be more pro- 
nounced in the case of a secular drift as the time baseline of the 
measurements increases, which is the case when ELODIE obser- 
vations are considered. The scatter also reflects the RV errors, 
which are higher for ELODIE than for the other instruments (see 
Table[T]i. Table|2]summarizes the mean difference and RMS scat- 
ter obtained when comparing the measurements from two instru- 
ments. In this process, an iterative clipping at the 3<x level is per- 
formed to remove the outliers. The offset column provides the 
velocity corrections applied to the measurements of CORALIE, 
NARVAL, ELODIE, and HARPS to translate them onto the SOPHIE 
scale. 

iBoisse et alJ ([2012) have also measured the RV offset be- 
tween SOPHIE and ELODIE and found a steeper slope (-425.6(5 - 
V) + 202.4) but their fit is limited to B - V < 0.75. Tableland 
Fig. [5] also present the comparison of the CORALIE and HARPS 
measurements as an illustration of the excellent agreement be- 
tween those two instruments. 

Two tables are provided for the catalogue, only available 
electronically. Table[3]gives the original individual RV measure- 
ments, and Table|4]gives for each star : 
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Fig. 5. Offsets between RV measurements from HARPS and 
CORALIE showing a remarkably low RMS scatter. 



the basic data from Hipparcos: HIP number, equatorial coor- 
dinates, V magnitude, B-V, spectral type; 
the catalogue(s) from which the star was selected 



iNidever et all (I2002I) . iFamaev et all (|2005[) . iNordstrom et al] 
(12004 or IAU standard (lUdrv et al.Ul999l) : 



the weighted mean radial velocity in the SOPHIE scale, RVs, 
with a weight w; applied to each individual velocity mea- 
surement RV, being w, = 1 / e' 2 and e' the quadratic sum of 
et and the offset error listed in Table |3j 
the internal error of RVs : / = Y* w i ^lEW, 
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Table 2. RV offsets between instruments derived for stars in 
common. 



Instruments 


N 

ly common 


N, 


Offset 






(a) 


(m s -1 ) 


CORALIE- SOPHIE 


164 


154 


-23 ± 2.4 


NARVAL-SOPHIE 


114 


108 


-26 ± 6.2 


ELODIE-SOPHIE 


277 


247 


(b) 


HARPS- SOPHIE 


36 


34 


17 ±5.0 


HARPS-CORALIE 


88 


84 


31 ± 1.5 



Notes, (a) N, is the number of stars effectively used after an iterative 
clipping at a 3cr level. 

(b) In the case of ELODIE-SOPHIE a linear regression was performed, 
giving an offset as a function of the colour index Offset= -259.0(6 - 
V) + 105.2 + 2.9 ms-'. 

- the RVs weighted standard deviation cr RVs defined as 

X>i 

4 Vs = (V L v 2 I>( RVi - ™*f 
(!>/) - I>; ; 

i I 

- the num ber of observations N; 

- the RVs uncertainty defined as the maximum of the standard 
error <x RVs / V/V and // V/V (Jasniewicz &Mavor[ fl9881: 

- the time baseline in days; 

- the mean Julian day of the N observations. 



Table 3. Excerpt of the catalogue table giving the 10214 original 
individual RV measurements and their error e (see Sect. [3]), the 
Julian day of the observation, the instrument, and mask used. 
The name of the instrument is abbreviated with its first letter. 



ID 


Instr. 


JD 


RV 


e 


mask 






-2400000 


km s 


km s 




HIP000296 


C 


52173 


10.932 


0.0070 


G2 


HIP000296 


c 


53262 


10.930 


0.0072 


G2 


HIP000296 


c 


54354 


10.952 


0.0064 


G2 


HIP000296 


c 


55167 


10.948 


0.0109 


G2 


HIP000407 


c 


51542 


12.346 


0.0093 


G2 


HIP000407 


c 


52649 


12.316 


0.0069 


G2 


HIP000407 


c 


53921 


12.328 


0.0075 


G2 


HIP000407 


c 


55167 


12.331 


0.0123 


G2 


HIP000420 


c 


51382 


34.861 


0.0082 


G2 


HIP000420 


c 


53218 


34.852 


0.0083 


G2 


HIP000420 


c 


55167 


34.868 


0.0130 


G2 


HIP000420 


c 


55758 


34.879 


0.0068 


G2 


HIP000466 


s 


55131 


-8.371 


0.0080 


K5 


HIP000466 


s 


55452 


-8.364 


0.0080 


K5 


HIP000556 


c 


55882 


0.109 


0.0110 


G2 


HIP000556 


c 


56221 


0.091 


0.0135 


G2 


HIP000556 


s 


55129 


0.102 


0.0086 


G2 



Figure [6] shows the histogram of the variations of the 1420 
candidates. We define the level of stability of a given star by 
3cr RVs . We find that 92.8% of the stars have a stability better 
than 300 m s" 1 , which is the threshold defined for the calibration 
of the RVS instrument. More than 1000 stars exhibit a stability 
better than 100 m s _1 (i.e. cr RVs < 33 m s~'). 

5. Comparison with other catalogues 

The catalogue of Nide ver et al.l (2002) contains 889 late-type 
stars followed during several years, with 782 of them exhibit- 
ing a velocity scatter less than 100 m s _1 . Some 336 stars of 




RV variation (m/s) 



Fig. 6. Distribution of RV variations of the candidate standard 
stars having at least two RV measurements separated by 100 
days or more. The variation is defined by 3cr RVs . A dashed line 
shows the 300 m s _1 stability threshold adopted for the calibra- 
tion of the RVS instrument. 



the catalo gue presented here have been selected from the stable 
stars of iNidever et al.l (|2002) and fulfil the Gaia-RVS require- 
ments. Our catalogue also has 3 55 stars in common with the 
recent list of IChubak et al.l d2012l) which provides RVs for 2046 
nearby FGKM stars. Comparisons with those two catalogues are 
shown in Fig. [7] as a function of B - V. There is general good 
agreement with small zero-point differences. The outliers reveal 
the intrinsic variable stars, as well as systematics due to method- 
ology. For B — V > 1.2 there is a cl ear systematic effect. We note 
that both INidever etal.1 (|2002|) and lChubak et all (1201 2|) use the 
same technique of chi-square minimization to determine the rel- 
ative shifts between a template spectrum and an observed spec- 
trum, but not in the same spectral range. They use an M-dwarf 
observed spectrum as template for cool stars and a solar spec- 
trum for the other stars. The t wo catalogues a re con structed on 
the same scale, as explained in Chuba k et alj d2012l) . We divide 
the comparison of their values to the present ones into two parts, 
according to their reference template. Results are presented in 
Tabled 

We get very similar zero-point differ ence in both cases , 
whic h is expected s ince the catalogues of INidever et al.l ([2002) 
and IChubak et al.l d2012l) are on the same scale. The very 
low RMS scatter of 33 m s obtained for the comparison to 
INidever et al.1 (120021) for FGK stars is remarkable for two rea- 
sons. First, since the median epoch of Nidever et al. observa- 
tions (1999.0) is earlier than ours (2006.9), it means that most 
of the stars in common have not varied much during that period. 
Second, it also shows that the precision of the present catalogue 
is not worse than 33 m s" 1 . The same level of scatter has al- 
ready been measured by INidever et al.l (120021) when they com- 
pared their catalogue to the best set of stand ard stars available 
at the time, assembled by the Geneva team dUdrv et all [l999) 
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Table 5. C omparison of RVs o btained with by iNidever et al.l 
(2002) and IChubak et alJ (120121) and those presented here, ac- 
cording to the templates they used. Differences, Arv = RVs - 
RVchubak/nidever (Fig. U} and RMS scatter in m s -1 are de- 
rived with an iterative clipping at a 3cr level. N* is the number of 
stars in common, N„ is the number of outliers removed. 



source 


stat 


Sun template 


M-dwarf template 




N. 


326 


29 


Chubak et al. 


N„ 
Arv 


15 

63 



-98 




RMS 


100 


182 




N, 


308 


28 


Nidever et al. 


N„ 
Arv 


22 
72 




-141 




RMS 


33 


178 



and based on CORAVEL and ELODIE. The ELODIE measurements 
that were used to build the Udry et al. standard star catalogue are 
among those retrieved from the arch ive for the p r esent c atalogue. 
We do not show the comparison to lUdrv et alJ (Tl999) since we 
have measurements in common. It is worth noticing that th ere is 
no correlation of the differences with lNidever et al.l d2002) with 
B - V in the FGK range, even though the stars of the present 
catalogue have been ana lysed with different m asks. 

Fo r RVs obtained bv lChubak et al.1 d2012l) and lNidever et alJ 
(2002) with an M-dwarf template, the mean difference and RMS 
are much larger. Figure |7]suggests that there might be a correla- 
tion of these differences with colour, but the B — V interval is too 
narrow to confi rm this cla i m. Su ch significant offset and RMS 
compared with lUdrv et ail ([1999) have already been found and 
discussed in bot h articles, so we con firm with a larger sample 
the assertion by Chuba k et al.l d2012l) that 'the M dwarf veloci- 
ties in general, from all surveys, remain uncertain at the level of 
200 m s^ 1 (RMS) and harbor uncertain zero points at the level 
of 150 ms- 1 '. 

The comparison with the catalogue of lLatham et al.l d2002l) is 
shown in Fig. [8] This catalogue contains 1359 single-lined stars 
on the Center for Astrophysics (CfA) scale. The difference is: 
Arv = RVs - RVlatham = 161 m s -1 for 66 stars in common 
after rejection of two outliers. This offset is quite significant but 
it agrees with the zero-point correction of 139 m s _1 that they 
determined from asteroids. Although the Latham et al. RVs are 
presented as having errors ranging between 0.5 to 1 .5 km s _1 , we 
measure a much lower RMS scatter, 164 m s _1 , which suggests 
that their errors are lower for the 65 common stars. 
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Fig. 8. Comparison of RVs obtained with by Latha m et al.l 
(2002) and those presented here. The offset and RMS values 
have been obtained after rejectiing two outliers. 
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Fig. 9. Difference in radial velocities of the Chubak et al. stan- 
dard stars for 63 stars in common. The offset and RMS values 
have been obtained after rejecting seven outliers at B - V > 1.2. 



6. Stable and variable stars 

IChubak et all d2012l) present a list of 131 FGKM standard stars 
having constant radial velocity over ten years. We have 63 stars 
in common with this list and the comparison is presented in 
Fig. [9] The average offset is Arv = RVs - RVchubak = 81 m s _1 , 
with an RMS scatter of 63 m s after rejection of seven outliers 
at B - V 1.2. This RMS scatte r is lower than in th e compari- 
son with the main catalogue of IChubak et al] d2012l) presented 
in Table[5] but does not reach the very low value of 33 m s -1 ob- 
tained in the comparison with lNidever et al. (2002). This means 
that either the Chubak et al. standard stars are less stable than the 
Nidever et al. ones, or that their errors are significantly higher 
than the Nidever et al. ones. 

The catalogue presented here is also a good source for very 
stable stars that are good candidates for becoming official RV- 



standards useful for many projects other than Gaia. Table 2]can 
easily be queried with criteria such as 

- at least 4 measurements available, 

- standard deviation on the SOPHIE scale cr RVs < 33 m s _1 (i.e. 
stability better than 100 m s~'), 

- time baseline of ten years at least. 

More than 300 stars verify these conditions and should 
be followed up in the future to become official RV-standards. 
Several of these stars are shown in Fig. [10] They were observed 
with different instruments and the agreement of the measure- 
ments demonstrates that the zero points have been properly cor- 
rected. 

Several stars showing a remarkable trend in their RV are 
shown in Fig. Q~T] The presence of a low-mass companion 
likely explains these drifts. The oscillation of HIP017960 is 
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due to a Jupiter like exoplanet already reported by Boiss e~eTafl 
(120121). HIP079248 is als o known to harbour two exoplanets 
(Gozdzie wski et al.L 120061) . A Keplerian orbit with a period of 
1161 + 45 days was fitted on HIP098714 measurements by 
ISivan et ail (1200 4). revealing the presence of a planetary com- 
panion. Its semi-amplitude of 95 ± 17 ms"' is still low enough 
to make this star suitable for the RVS calibrations. Interestingly, 
HIP041844 and HIP044089 are part of the Nidever et al. stable 
stars. The drift of these two stars is quite slow, which might ex- 
plain why Nidever et al. have not been able to detect it within the 
time range of their observations. This demonstrates that it is very 
important to have observations over a sufficient time base to de- 
tect such variations. However, the number of binary stars remain- 
ing in our samp le is expected to be very low. The selection crite- 
ria described in ICrifo etal.l (l2010h and recalled in Sect. |2]elimi- 
nated all known double stars and variable stars (photometrically 
and spectroscopically), as well as those with fainter neighbours. 
HIPPARCOS efficiently detected double stars with separations 
of 0.1" to 10", and magnitude differences lower than 4 magni- 
tudes (see Fig. 3.2.106 available in the HIPPARCOS Catalogusfl 
volume 1 ). Su ch stars are labell ed in HIPPARCOS and were 
not selected bv lCrifoetail(l2010l) . This guarantees that there is 
no spectroscopic binaries with short-term variations in our sam- 
ple, because the candidates are very nearby stars, all closer than 
150pc except some giants. Only a few wide systems with a faint 
secondary (M star or white dwarf) and a rather long period may 
therefore remain in the sample. For those, the second round of 
observations during the Gaia mission will be decisive in con- 
straining their long-term behaviour. With new observations to be 
done after 2014, we will extend the time baselines by several 
years (current median epoch of our observations is 2006.9), and 
we will have at least three measurements per star (67% of the 
stars have already 3 measurements or more), well distributed in 
time and with a much higher precision than the threshold of 300 
m s -1 . This will guarantee that we cannot miss any long period 
binary with large amplitude that would not be suitable for the 
RVS calibrations. 

It is interesting to investigate whether the most stable stars, 
with <t RVs < 30 m s , or on the contrary the less stable stars, 
with <x RVs > 100 m s , display in peculiar locations of the HR 
diagram. Figure [12] shows the HR diagram of the 1420 candi- 
dates with a colour code corresponding to their stability level 
3<r RVs . The median <x RVs and the fraction of stable and non- 
stable stars are evaluated into four parts of HR diagram, as in- 
dicated in Fig. [12]: main sequence (My > 4.5, B - V > 0.65), 
turn-off (B-V < 0.65), giants (M v < 2.5, B-V > 0.85), and sub- 
giants (the rest of the sample). Results of the counts are given 
in Table [6] They clearly show that giants are less stable on av- 
erage than dwarfs and subgiants, confirming previous findings 
(Dollinger et all 120051) . Not only is the presence of a stellar or 
substellar companion orbiting around the giant star suspected to 
cause the RV variations, but also the possible rotational modula- 
tion of surface inhomogeneities or pulsations. 

7. Next steps 

There are many physical mechanisms that af fect spectroscopic 
radial velocities (Lindeg ren & Dravinsl 120031) . For instance, the 
convective shifts caused by motions in stellar atmospheres de- 
pend on stellar lines and on the temperature and gravity. They 
can reach 3 km s _1 for a Can line in an F dwarf, and -0.4 km s _1 

2 on-line at 

www . r ssd . esa . int/index . php?pro j ect=HIPPARC0S&page=0verview 
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Fig. 12. HR diagram of the 1420 candidates with a colour code 
corresponding to their stability level 3cr RVs . The less stable stars, 
with 3<x RVs > 300 m s _1 , are indicated as black points. Those 
are not suitable for calibrating the RVS instrument. The HR dia- 
gram is divided into four parts corresponding roughly to turn-off, 
main sequence, subgiants, and giants where median RV scatter, 
and the fractions of stable and variable stars have been evaluated 
(Table©. 



Table 6. Percentage of stable (cr RVs < 30 m s ) and variable 
(c R v s > 100 m s~') stars among our 1420 selected candidates, 
in four parts of the HR diagram. 



HR part 


N 


median tr RVs 


stable 


variable 






m s 


% 


% 


main sequence 


397 


17.0 


76 


9 


tum-off 


662 


18.3 


74 


6 


subgiants 


199 


17.7 


75 


5 


giants 


154 


34.6 


46 


11 



for an Fei line in a K giant (IChiavassa et all 1201 ll) . Other as- 
trophysical processes may affect the spectroscopic RV of a star, 
such as its rotation and activity, the granulation, or the presence 
of low-mass companions. 

One of the objectives of our ongoing project will be to pro- 
vide an accurate kinematical RV for each of the 1420 stars in 
the sample. This will be done by deriving the gravitational red- 
shift from atmospheric parameters and the convective shift cor- 
rections from three-dimensional hydrodynamical model atmo- 
spheres. Another step will be to determine the zero point of the 
catalogue by comparing the RV measurements of a selection of 
asteroids with their kinematical velocities as derived from celes- 
tial mechanics. Finally, we will continue the observations of the 
1420 candidate stars during the Gaia mission in order to validate 
the long-term stability of their RV. All those steps are essential 
for calibrating the measurements that will be performed with the 
Gaia-RVS. 

8. Conclusion 

We have presented the pre-launch version of the catalogue of 
RV standard stars for Gaia, assembled thanks to a long-term ob- 
serving programme started in 2006 on several spectrographs and 
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with archived measureme nts. The precis i on of this catalogue is 
at the same level as that of Ni dever et al.l (120021) . ~ 33 m s _1 , for 
FGK stars (0.5 < B - V < 1.2), showing that the RV measure- 
ments from the different instruments have been properly com- 
bined. The vast majority of the 1420 selected candidates are 
found to be stable at the 300 m s level, which makes them 
suitable for calibrating the RVS instrument. Their long-term sta- 
bility will have to be confirmed with new ground-based observa- 
tions during the Gaia operations. 

The stable stars presented here have characteristics that make 
them very useful as standards for many projects other than Gaia. 
They have a good coverage of important parameters, such as sky 
distribution, apparent magnitude, spectral type, and luminosity 
class. They can be used to intercompare intruments and calibra- 
tions of zero point of spectroscopic observations. The stars and 
measurements of this list and those of all other studies on RV- 
calibrators will soon be discussed within the task group 'radial- 
velocity standard stars' of the IAU Commission 30. 
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Fig. 4. Offsets between RV measurements from SOPHIE and the other instruments as provided in Table [2] Error bars when available, 
are the quadratic sum of the two standard deviations obtained for a given star observed several times with the two instruments being 
considered. 



Table 4. Excerpt of the catalogue table giving the mean RV on the SOPHIE scale of the 1420 standard star cand idates (see text for th e 
de scription of the colum ns). Abbreviations are used fo r the catalogue(s) from which the star was selected: ni= lNidever et al. (2002), 
fa= lFamaev et all (120051. nor- Nordstr om et al.l (12004 in column 6. 



ID a, <S Vmag B-V ST cat RVs I o"rv s N uncertainty span mean JD 
bus' 1 kms' 1 kins'' kms' 1 days -2400000 



HIP000296 


00:03:41.30 -28:23:45.1 


8.24 


0.780 


G8V 


nor 


10.963 


0.0078 


0.0117 


4 


0.0059 


2994 


53739 


HIP000407 


00:04:58.63 -70:12:43.9 


8.13 


0.710 


G5V 


nor 


12.351 


0.0086 


0.0130 


4 


0.0065 


3625 


53319 


HIP000420 


00:05:07.60 -52:09:05.1 


7.53 


0.577 


GOV 


nor 


34.889 


0.0086 


0.0129 


4 


0.0065 


4376 


53881 


HIP000466 


00:05:34.91 +53:10:18.1 


7.23 


1.174 


K0 


fa 


-8.367 


0.0080 


0.0049 


2 


0.0057 


321 


55291 


HIP000556 


00:06:46.94 -04:20:59.3 


8.21 


0.574 


F8 


nor 


0.113 


0.0104 


0.0165 


3 


0.0095 


1092 


55744 


HIP000616 


00:07:32.34 -23:49:08.2 


8.70 


0.798 


K0V 


ni 


-42.994 


0.0086 


0.0049 


4 


0.0043 


3380 


53629 


HIP000624 


00:07:37.41 -45:07:09.5 


8.19 


0.584 


F8/G0V 


nor 


14.802 


0.0117 


0.0134 


2 


0.0095 


725 


55145 


HIP000699 


00:08:41.02 +36:37:38.7 


6.21 


0.504 


F8IV 


ninor 


-15.116 


0.0043 


0.0119 


16 


0.0030 


5160 


54912 
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Fig. 7. Difference between our radial velocities and those from iNidever et alJ J2002) (left panel) and Chubak et al] (I2012I) (right 
panel) as a function of the B-V colour index. Red dots indicate the stars they have analysed with an M template, while for the other 
ones they used the solar spectrum as reference. The offset and RMS values are given for FGK stars and for M stars in Table [5] 
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Fig. 10. RV measurements for some of the long-term stable stars of the catalogue. The RV axis is centred on RVs and spans 300 
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Fig. 11. Trend in RV measurements for some stars of the catalogue. Symbols as in Fig.fTUl 
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